Quantum metamaterials: Electromagnetic waves in a Josephson 

qubit line 

A.L. Rakhmanov 1,2 , A.M. Zagoskin 1 ' 3 ' 4 , Sergey Savel'ev 1,3 , and Franco Nori 1,5 
(1) Frontier Research System, The Institute of Physical and Chemical Research (RIKEN), 

Wako-shi, Saitama, 351-0198, Japan 
(2) Institute for Theoretical and Applied Electrodynamics RAS, 125412 Moscow, Russia 
(3) Department of Physics, Loughborough University, 
Loughborough LE11 3TU, United Kingdom 
(4) Physics and Astronomy Dept., The University of British Columbia, 
Vancouver, B.C., V6T 1Z1, Canada and 
(5) Department of Physics, Center for Theoretical Physics, 
Applied Physics Program, Center for the Study of Complex Systems, 
University of Michigan, Ann Arbor, MI 48109-1040, USA 
(Dated: February 1, 2008) 

Abstract 

We consider the propagation of a classical electromagnetic wave through a transmission line, 
formed by identical superconducting charge qubits inside a superconducting resonator. Since the 
qubits can be in a coherent superposition of quantum states, we show that such a system demon- 
strates interesting new effects, such as a "breathing" photonic crystal with an oscillating bandgap, 
and a "quantum Archimedean screw" that transports, at an arbitrary controlled velocity, Joseph- 
son plasma waves through the transmission line. The key ingredient of these effects is that the 
optical properties of the Josephson transmission line are controlled by the quantum coherent state 
of the qubits. 
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I. INTRODUCTION 



The development of superconducting electronics now allows the observation of quantum 
behavior, such as the coherent superposition of different macroscopic states in meso- and 
macroscopic devices^. While efforts in this field are now mainly directed at the development 
of superconducting quantum bits as elements of eventual quantum computers, there are other 
interesting possibilities opened by the existence of such large, controllable, quantum coherent 
circuit elements. In particular, analogies with cavity QED have led to interesting theoretical 
and experimental results&^&Z^. There, qubits play the role of artificial atoms, and high- 
quality superconducting resonance circuits mimic optical cavities. Differences include that 
instead of having a stream of identical atoms moving through the cavity, the state of a single 
qubit, permanently coupled to the resonator, could be periodically changed. In Ref. I] results 
of an earlier experiment^ were considered from the point of view of lasing in such a system. 

Another recent surge of interest in the electrodynamics of Josephson Junction arrays is 
related to THz electromagnetic waves propagating in such systems. THz waves are important 
for applications, but are hardly controllable for both optical and electronic devices. Thus, 
Josephson structures might be of potential importance for miniature THz generators, filters, 
detectors, and wave guides- 11 . 

In this paper, instead of considering a single or a few qubits, we investigate the behaviour 
of an infinite chain of identical qubits inside a resonator, from the point of view of THz or 
sub-THz electromagnetic wave propagation in such a quantum medium. We show that by 
placing the qubits in a quantum superposition state, some interesting possibilities can be 
realized, including "breathing" photonic crystals and an "Archimedean screw" transport of 
classical electromagnetic modes. 

For the lack of a better term, we call such qubit structures, considered from the point 
of view of macroscopic propagation of electromagnetic field, quantum metamaterials. This 
because (classical) metamaterials allow additional ways to control the propagation of elec- 
tromagnetic fields, not available to standard materials. (Alternative approaches to super- 
conducting metamaterials were investigated in Refs. QQQQQQQ 

.) Similarly, our 

proposed quantum metamaterials allow additional ways of controlling the propagation of 
electromagnetic waves, not possible with normal classical structures. Indeed, the coherent 
quantum dynamics of qubits determines the THz "optical" properties in the system. 
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II. MODEL 



As a model, we choose a set of identical charge qubits placed at equal intervals, /, between 
two bulk superconductors separated by a distance D (Fig. [T]). Each qubit is a small super- 
conducting island connected to each superconducting bank by a Josephson junction. The 
superconducting phase on the nth island is <p n . When treated quantum mechanically, such 
an island indeed consitutes a qubit, if its total capacitance is small enough^. The magnetic 
field H is applied normal to the structure (in the y direction) and the vector potential A has 
only a z-component. We denote by A zn the vector-potential between the n-th and (n + l)-th 
qubits. 

The structure of Fig. [T] is a ID waveguide with the energy per unit length 
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Here the dot denotes d/dt, the Josephson energy, Josephson frequency, critical current and 
junction capacity are respectively Ej = Q I c /2tic, uj = eI c /hC, I c and C\ $o = h/2e 
is the flux quantum. We took into account that, in the presence of the vector potential, 
the superconducting phase differences across the junctions of the nth qubit, ±yj n , acquire a 
gauge term, a n = 2nDA zn /^> . Introducing the dimensionless units E = S/Ej, and t — > u t, 
we rewrite Eq. (JTJ) as 



E = (p n + a n — 2 cos a n cos (f n + (3 (a n+i - a n ) 



(2) 



where 



E 



EM 



(3) 



SnDlEj \2n J Ej 
characterizes the ratio (E EM /Ej) of electromagnetic and Josephson energies. 

In this paper the electromagnetic (EM) field is treated as a classical wave. We also assume 
that its amplitude is small, a n < 1. This means that the magnetic flux per unit cell area 
H y D x I is much smaller than <3> - Under such assumptions, the Hamiltonian for a single 
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FIG. 1: (Color online.) Geometry of the system. Identical charge qubits are placed at equal 
intervals I between bulk superconductors separated by a distance D. 



We restrict the states of each qubit to either its ground state |0), with energy E , or excited 
state |1) with energy E\. This is justified due to the nonlinearity of the Josephson potential. 
We are not concerned here with decoherence effects in the qubits, concentrating exclusively 
on their interaction with the electromagnetic wave in the system. This idealization is justified 
as long as the decoherence time exceeds the wave propagation time across a significant 
number of unit cell periods. This is not unrealistic, given the high quality of superconducting 
resonator-qubit devices already achieved by recent experiment (e.g., in Ref. If] the quality 
factor Q exceeded Q = 10 4 at 6 GHz). 

In the absence of an EM field, the wave function of the system, \I/ ra , is a sum 
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where C£° are constants and e is the dimensionless excitation energy, 

_ E x — E 

In the presence of an EM field, the coefficients C% become time-dependent, and, as it follows 
from Eq. (jlj), these obey the relations 19 

i ^r = a " £ v km(t) cm (6) 

m=l,2 

with the initial conditions C£(t = 0) = C^°. Here 

CW = (A;|cos^ n |m) 

are matrix elements of the interacting field-qubit interaction from (jlj), calculated in the 
Heisenberg basis 

{ |0) exp (iet/2), |1) exp (-iet/2)} . 

Later we will also use the time-independent matrix elements V£ m = (k\ cos (p n \m) in the 
basis {|0), |1)}. 

Varying the energy (T2]), we obtain the equation for the electromagnetic field in the linear 
approximation: 

a n - f3 2 (a n+ x + a n _i - 2a n ) + a n ($ n \ cosip n \^ n ) = 0. (7) 

The set of Eqs. ()5]), ([6]), and ([7]) should be supplied by appropriate initial and boundary 
conditions. By controlling the qubits in Eq. (Q, we propose to change the transmission and 
reflection of EM waves described by Eq. (171). 

We are interested in the case when the wavelength is large compared to the size of the unit 
cell. Therefore the qubit line can be treated as a continuous ID medium with n-l replaced by 
x. The difference equation d7J for a n (t) and ^ n (t) is thus replaced by a differential equation 
for a(x, t) and ^(x, t) 

& ~ + V * a = °> V ° = (*( x )\ cos V (x)\V(x)) (8) 

Within the perturbation theory approach, we present the electromagnetic wave as a sum 
of the larger incident wave, cto, and a smaller scattered wave a±. A quantum state of the 
system is described by the wave function 

V(x,t) = C (x,t)|0)e ie * /2 + C 1 (x,t)|l)e- i£ * /2 . (9) 
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In the unperturbed state the coefficients in this equation are Cj = Cf(x). We present the 
coefficients Ci(x,t) as a sum of the unperturbed solution Cf(x) and a small perturbation 
Cftx,t), Ci(x,t) = C°(x) + Cl(x,t), with \C}\ < 1. Using Eq. ©, we derive 

iC\ = J* dt'al (V 00 C ° + V 01 C?e- fc *') 

itfj = jf dt'al (y n C? + VTo C °e l£ '') , (10) 

where = (i| cos 99 1 A;) are calculated using the unperturbed wave functions, V* means 
complex conjugate of V, and V* Q = V 01 . 

For the unperturbed EM wave ao, we obtain from Eq. (jSj) 

«o-/5 2 ^ + K.«o = 0. (11) 

Here, (Vo) 1//2 plays the role of the Josephson plasma frequency, which is now controlled by 
the quantum state and quantum dynamics of the qubits. For the matrix element Vq we can 
derive the following expression 

Vq = |C °| 2 V 0Q + \C° X \ 2 V u + CgCf e* 1 V w + h.c. (12) 

For simplicity, we assume that «o is a standing wave, «o = Acos(ujt) cos[k(cu)x]. 



III. ELECTROMAGNETIC WAVE PROPAGATION THROUGH A UNIFORM 
QUBIT LINE 

A. Qubits initially in the ground state |0). 

If all the qubits are in the ground state |0), then initially Cq = 1 and = 0. In this 
case Vq = Vqq and the wave vector is 

k(u) = ^Vu 2 -V 00 . (13) 
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Thus, the wave can propagate if its frequency exceeds (Voo) 1//2 , which can be interpreted as 
the "ground state" plasma frequency of the medium. From Eq. (TlOl) we obtain 

C£(x,t) iA 2 cos 2 (kx) f sin(2u;t)| 



f t , sin(2^t) j 



Voo 2 
Cl(x,t) __ A 2 cos 2 (kx) J e i£t - 1 
Voi : 2 \ e 

e + e i£ * [2iw sin(2cut) - £ cos(2u;t)] 



(14) 



+ 



Au 2 - e 2 



The initial disturbance of the wave function produces a disturbance a± in the propagating 
wave. For this perturbation, using Eq. (jSJ), we derive 



Si - /? 2 ^ + Hoai + AVoa = 0, (15) 
AV being the perturbation of the field-qubit coupling. By means of Eqs. (jl~4|) we find 

\V 01 \ 2 A 2 cos 2 {kx) 

(16) 




2(2u; 2 - e 2 ) cos(et) + e 2 cos(2ujt) 



e e (4oo 2 - e 2 ) 

We see that the electromagnetic wave is in resonance with the qubit line if its frequency is 
half the inter-level distance, u = e/2. This is due to the term proportional to a 2 in the 
Hamiltonian (j4j). Near the resonance, the condition \C}\ <C 1 is no longer valid and the 
usual perturbation approach fails. 



B. Qubits initially in the excited state |1). 

If all qubits are initially in the excited state |1), the solution is found in complete analogy 
to the previous case. As it can be readily seen, we should only exchange «-> 1 and e «-> —s 
in Eqs. (TI5])- (fIBl) 

k(u) = ^y/^-Vu, (17) 



C\(x,t) iA 2 cos 2 (kx) ( sin(2c;t)l 
V n ~ 2 [ t+ 2^ J 

Cl(x,t) A 2 cos 2 (kx) \ l-e~ iet 
Wo ' 2 | i 

-e + e" i£ * [2zcj sin(2cjt) + e cos(2a;t)] 
+ Au 2 - e 2 



(18) 



and for the electromagnetic wave we obtain 



on — ft 




+ V n ai + AVia = 0. 



(19) 



dx 2 



In Eq. (fT9|) 



AV^t) 



AV (t) 



where AVo(i) is given by Eq. ( fl6l) . 

The excited qubit line is an active medium, and one should expect a resonance pumping 
of the electromagnetic wave as it propagates along. However, AVi —>■ at 2u — > e (and 
AVo — > at 2u — > £ as well). This "paradox" only reflects the limitations of the first order 
perturbation approximation, where |C° + Cf\ 2 = | C ? | 2 . In other words, to first order, the 
qubit energy does not change. To describe the pumping effect, we must take into account 
the higher order terms, which is beyond the scope of our current investigation. 

Finally, let Cq = C\ = 1/2. In this case all qubits 'rotate ' between the ground and excited 
state. The matrix element in Eq. (jSJ) is now 



Let us now assume, for simplicity, that the frequency of the electromagnetic wave is high, 
uj e. Then its wave vector is a slowly oscillating function 



then the qubit line will alternate between transparent and reflecting state with a frequency 
e, as the wave vector k(t) switches between real and imaginary values. In addition, the qubit 
line produces electromagnetic waves with frequencies e and u ± e. 

IV. ELECTROMAGNETIC WAVE AT RESONANCE 

Near the resonance, 2uj — e = Auj <C oo, we can use a resonant perturbation approach^, 
instead of a first-order perturbation approximation. In doing so, we drop all the terms in 



Vo(t) = 7 [Voo + V u + 2V m cos(et)] . 



(20) 





If the wave frequency u is close to the threshold value, 



yVoo + Vu/20, 
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Eqs. (E} except the resonant ones. As a result, we have 



iC = -Vle tAuJt d, 
iC l = -Q*e~ lAut C , 



where 

The solution of this system is 

e iAwt/2 



C (t) 



n* 



Q(x) = A 2 cos 2 [k{u)x]V Q1 /4. 



^ A -y) e8Ai -^( A+ y) e " !Ai 



(22) 



C x {t) = e - i ^ t ' 2 ( 9l e iAt + g 2 e- iAt ) 



(23) 



where 



A = y fi 2 + (Alo) /4 , 

and gi, g 2 are constants. The coefficients Cj satisfy the normalization condition | C 1 2 + 
|Ci| 2 = 1. If at t = the system was in the ground state, we obtain 



C {t) 
Ci(t) 



JAu>t/2 



m\ 



. , . iAco . , . . 
cos(At) - — — sm(At) 



-iAut/2 



A 



2A 



sin (At). 



(24) 



At the resonance Au;=0, 



#=exp (iet/2) cos(|0|t) |0>— i exp (-iet/2) sin(|fi|t)|l), 

and each qubit periodically oscillates between its ground and excited states. The frequency 
of these transitions varies with the qubit position, since Q = Q(x). It can be considered as a 
spatially-dependent quantum beat frequency. If the system was initially in its excited state, 
then, we have 



C (t) 
Ci(t) 



iAujt/2 



A 

-iAwt/2 



sin (At), 



, , , iAu . . , . 
cos(At) - — — sm(At) 



(25) 



and at Au;=0, 



^=^exp (iet/2) sin(|fi|t)|0)+exp (-iet/2) cos(|f2|t)|l). 
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Using these expressions for the wave functions, we find the matrix element Vq in Eq. (fl2l) . 
Assuming for simplicity that Vn~ are real and Auj 2 <C |f2| 2 , we obtain 

Vo(t) = g ± 2 COS ( 2 I°IO 

- V r isin(2|O|t)sin(£* + Aw), (26) 

where the signs '+' and ' — ' correspond to the initial ground and excited state for qubits, 
respectively. If the transitions between the ground and excited states are suppressed, | Voi | <C 
| Vo — Vii|, then, in resonance, the electromagnetic wave a has a time-dependent wave vector, 

k(u,t) = y/u*-V Q {t), 

where Vo(t) varies between Voo and Vu. These results are valid if the wave-qubit interaction 
does not distort it too much: \a\\ <C \o$\. In any case the condition 

\n\ ~ A 2 \V 01 \ < u ~ e (27) 

must be fulfilled. 



V. QUANTUM METAMATERIAL PHOTONIC CRYSTAL 

In analogy with photonic crystals^, the interaction of the electromagnetic wave with 
qubits can produce a frequency gap in the spectrum of the propagating wave, if the qubit 
states are periodically modulated in space. For example, suppose the qubits are in the (7) 
or 1 5) state with a spatial period 2L. The wave obeys the equation 

a - (3 2 a xx + V~ r/ a = (28) 

or 

a - (3 2 a xx + V 55 a = 0. (29) 

The states I7) or \5) can be either stationary (eigenstates of the qubit Hamiltonian) , or their 
superpositions. In the latter case, the photonic crystal discussion makes sense only if the 
quantum beat frequency is small compared to the frequency of the propagating wave, that 
is, e 2 < I Voo I and |Vn|. 

Following the usual band-theory approach for electrons in a crystal lattice, we seek the 
solution of Eq. (|28|29p in the form of a Bloch wave a(t, x) = u(x, k) exp(ikx — iut), where 
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FIG. 2: (Color online.) Photonic crystal spectrum obtained from a qubit transmission line. 
u(k) (red) shown was calculated for a periodic array of qubit states (ground and excited), 
ratio Vqo/Vu = 5. 



The 
The 



u(x, k) is a periodic function of x with the period 2L, and the dimensionless wave vector k is 
in the first Brillouin zone, —n/L < k < n/L. Consider the jth elementary cell of our periodic 
structure: for Xj < x < Xj + L all the qubits are in state (7), and for Xj + L < x < Xj + 2L 
in the state \S). In both regions, the solution a(t,x) of Eq. ( 12811291) is a sum of exponential 
terms multiplied by constants Cj. Using the continuity of a and da/dx at the boundaries 
of different regions and the periodicity of the Bloch functions u(x,k), we obtain a set of 
homogeneous linear equations for Cj. The nontrivial solution of these equations exists only 
if the determinant of the set of equations is zero. Then, after straightforward algebra, we 
obtain the dispersion equation for the frequency u(k) in the form 



cos(k 7 L) cos(kjL) 




sin(/t 7 L) sin^L) 



(30) 
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FIG. 3: Breathing photonic crystal: contour curves of the wave vector k as a function of u and 
et in the situation described by Eq. (f34j) . The parameters used here are Vbo = Voi = 1, V\\ = 2, 
(3 = 0.5, L = 2. The time-dependent gaps in the spectrum are clearly seen. 



where 



77 



Vss 



(31) 



7 ^2 ' a p2 ■ 

This equation predicts the spectrum u(k) with gaps if the difference between k 7 and k$ is 
large enough; that is, \kL — Kg\ ^ 1, or 



IV- 
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V S5 \ >(3 2 . 



(32) 



Thus, in order to form a photonic crystal in the qubit line, the Josephson energy Ej must 
be large compared to the magnetic energy or, according to Eq. (J3j), 

2 



Ej > — 



$ 



(33) 



8tt I, 

The characteristic dependence of u(k) is shown in Fig. [2J Note that the gap value for the 
first zone is of the order of unity if the condition Eq. (|32l) is valid. 

The gap depends on the quantum state of the qubits, making this a quantum photonic 
crystal. Changing the microscopic quantum state of the qubits changes the macroscopic 
electromagnetic response of the system. 
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A more interesting situation arises if one or both of the qubit states are not the eigenstates 
|0) or |1), and therefore this system would exhibit quantum beats between the two. Let, e.g., 
|7) = |0) and 

\S) = {\0)e i£t + \l)e~ l£t } /2. 

Then V^ 7 = Voo and 

V ss (t) = [V 00 + V n + 2V 01 cos(et)]/4. 

In this case, the photonic crystal arises if any of the matrix elements is of the order of unity. 
The frequency gap is modulated by the value of Voi/2 with the period At = 2n/e. If Vqi ~ 1, 
then the modulation is significant. If 

| 7 ) = {\0)e i£t - \l)e- i£t ] /2 

and 

\5) = {|0)e ie ' + |l)e- i£ '}/2, 

then 

V 77 {t) = [V 0Q + V u - 2V i cos(e*)]/4, 

(34) 

V ss (t) = [Voo + V n + 2V 01 cos(et)]/4. 

In this case the photonic crystal appears if Vqi ~ 1. The gap is strongly modulated in time 
from zero, at t — (2n + l)7r/2e, to its maximum value, when t = nir/e (here n is an integer). 
Thus we obtain an interesting uniformly "breathing" photonic crystal, shown in Figj3l 

VI. QUANTUM 'ARCHIMEDEAN SCREW 

Let an external source produce a slow "control wave" propagating along the qubit line. 
That is, let the coefficients in the wave function be 

C (x, t) = sin (u t — k x), Ci(x,t) = cos (u t — k x), (35) 

where uq and ho are the frequency and the wave vector of the control wave, respectively, and 
both these values are small compared with the electromagnetic wave frequency and wave 
vector. Such a wave can be produced, e.g., by applying to the qubits a RF signal with a 
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FIG. 4: (Color online.) Quantum 'Archimedean screw': (top) schematic diagram of how the wave 
amplitude A(x,t)/Ao varies with position at a given moment in time; (bottom) the dependence 
of the wave amplitude on time and position; Vqo + V\\ = 2, Voo — V\\ = 1, (3 = 0.15, ko = 0.1, 
uj = 10- 3 , u = 1.23. 

position-dependent phase, and inducing Rabi oscillations between their ground and excited 
states. 

The matrix element Vq in the wave equations then takes the form 

v o{ x > t) = cos [2 (u t - k x)\ 

+ Voi sin [2 (co t — kox)} cos et, (36) 

where we now assume that Vqi is real. To simplify the problem, we now assume that 
| Voi | -C |Voo±Vii| and the last term, which describes the qubit relaxation, can be neglected. 
Thus, Vq = Vq(£), where £ = u t — kox. 
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The function Vo(£) varies slowly and we can seek the solution of Eq. (jSj) in the form 

a(t,0 = A(0 sm[u t- V (O}, (37) 

where the wave amplitude A(£) varies slowly and we can neglect terms with A"(£) (here 
/' = df/dt;). Substituting Eq. f|3T|) into Eq. (jSj) and separating terms with sin[a;ot — 
and cosfc^o^ — v(£)]i we derive two equations for the wave amplitude A and the phase r\ in 
the form 

2A> _ {(3 2 k 2 - u 2 ) r," 



(77 ) + 2ry - 



+ 



0. 



These equations are valid if 



A" 



A 



« (1/) 



(38) 



(39) 



A" 



(w - a; V) z 



(^W)' 



Integrating the first of these equations we get 

2 In A = - In (7/ + 



+ const. 



PS-W (40) 

The second equation is a quadratic equation for rf . Choosing its positive root and substi 
tuting it into Eq. (140]) we get 



A 



2\2 



mi - ui) 



1/4 



(41) 



(3 2 k 2 u 2 - {(3 2 k 2 - cu 2 ) V 
where A is a constant defined by the boundary conditions. The dependence A(x,t)/A is 
shown in Fig. H] for different values of parameters chosen within the validity range, Eq. (1391 ). 
It is easy to see that the wave amplitude A(t) achieves its maximum if the denominator in 
Eq. ()4ip becomes small at some moments in time. The value fi 2 k\Lo 2 — {(3 2 k 2 — tu 2 ) Vq should 
be positive. Assuming that Vqo > Vn, we can write the last condition in the form 



2 ^ 2 

ijJ > UJZ 



Vn[l- 



P 2 k 2 



(42) 



The variation of A(x, t) /A is maximum if u is close to u c ; however, the conditions in Eq. (I3"9"|) 
should be fulfilled. 

Returning to Fig. HI we see that the maxima of fast-oscillating electromagnetic field are 
transferred through the system at a pace, set by the much slower control frequency, u , 
reminding of the Archimedean screw, or a meat grinder. 
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VII. CONCLUSIONS 



Here we show that, for a classical electromagnetic wave, a line of qubits inside a super- 
conducting cavity plays the role of a ID transmission line with interesting characteristics 
and rich new physics. In particular, the quantum superposision of qubit states produces a 
"breathing" state with transparency changing with the quantum beat frequency of a single 
qubit. More interestingly, a periodic arrangement of qubit states yields a quantum photonic 
crystal, which can also be put into a breathing mode. A time-domain control of the qubits 
allows to realize an "Archimedean screw" state, where the incident electromagnetic wave 
is periodically modulated, and the regions of its maximum amplitude are carried along the 
qubit line with a desired speed. The investigation of the action of this system as an active 
medium requires a special treatment beyond the lowest-order perturbation theory, and will 
be the subject of future research. 
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